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The violation of mirror symmetry in the weak force provides a powerful tool to study 
the internal structure of the proton. Experimental results have been obtained that address 
the role of strange quarks in generating nuclear magnetism. The measurement reported 
here provides an unambiguous constraint on strange quark contributions to the proton's 
magnetic moment through the electron-proton weak interaction. We also report evidence 
for the existence of a parity-violating electromagnetic effect known as the anapole moment 
of the proton. The proton's anapole moment is not yet well understood theoretically, but it 
could have important implications for precision weak interaction studies in atomic systems 
such as cesium. 
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In 1933 the German physicist Otto Stern discovered that the magnetism of the proton was 
anomalously large, a factor of three larger than expected from the basic theory of quantum 
mechanics. This experiment turned out to be the first glimpse of the internal structure 
of the constituents of the atomic nucleus, and a tantalizing hint at the existence of quarks. 
Widespread applications of the proton's magnetic properties, such as the magnetic resonance 
imaging (MRI) techniques used in biology and medicine, have been developed despite a lack 
of fundamental understanding of the basic dynamics that generates this magnetism. After 
the key discovery of internal structure in the proton in a high energy electron scattering 
experiment at the Stanford Linear Accelerator Center in the late 1960's (1), the theory 
of Quantum Chromodynamics (QCD), which describes the interaction between quarks and 
the gluons that bind the quarks into the atomic nuclei observed in the periodic table, was 
developed. Despite almost 30 years of intense theoretical effort, QCD has been unable to 
produce numerical predictions for the basic properties of nucleons such as their degree of 
magnetism. 

In 1988, Kaplan and Manohar (2) proposed that the study of the weak magnetic force 
(analogous to the usual magnetic force associated with electromagnetism) would allow a 
separation of the proton's magnetism into the three contributing fiavors of quarks: up, 
down and strange. Measuring the contribution from strange quark-antiquark pairs is of 
special interest because it relates directly to the "sea" of virtual quark-antiquark pairs in 
the proton, a phenomenon predicted by QCD. In 1989, it was noted that the weak magnetic 
force could be isolated using its unique property of lack of mirror symmetry, or parity 
violation (3). The basic idea was to study the preference for the proton's interaction with 
electrons that have spin counterclockwise, over those with clockwise spin, relative to their 
direction of travel. We report data obtained using this method that, when combined with our 
previously published results (4,5), allow the first unambiguous determination of the proton's 
weak magnetism. We also report a measurement of a parity-violating, time-reversal-even 
electromagnetic contribution to proton structure, referred to in the literature as the proton's 
anapole moment. 

The weak magnetic moment of the proton, designated as iiz, is associated with the 
interaction due to exchange of the weak neutral Z particle. The weak and electromagnetic 
forces are related through the mixing angle dw, a fundamental parameter of the standard 
model of electroweak interactions that relates the Z boson to its charged counterparts, the 
(we use sin^ ^^=0.23117 ± 0.00016 (6) in our analysis). As a result, one can write nz 

as 

fJ'Z = (a*p - iin)-4: sin^ 9w jJ-p - jJ-s (1) 

where [Xp[n) are the usual magnetic moments of the proton (neutron) (2.97/ijv and — 1.91/iAr, 
respectively, where jiN = eh/2m.p = 3.152451238(24) x 10"^^ McV-T (6)). The individual 
contributions of up, down and strange quark-antiquark pairs to the proton's magnetic mo- 
ment, iiu, IJ'd, and Us, are defined by the relation Hp = l/^u ~ |A*d ~ |a*s • The neutron's 
magnetic moment is similarly constructed, interchanging only and fXd. Thus, a measure- 
ment of fiz combined with the known /Xp and provides the third observable required 
to uniquely determine the u, d and s quark pieces of the proton's magnetic moment. The 
study of parity-violating electron-nuclcon (e-A^) scattering enables a determination of the 
weak magnetic form factor Gf^{Q'^), and as a result the equivalent strange piece G%j{Q^), 
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the momentum-dependent counterpart of /i^, where is the relativistic four-momentum 
transfered to the proton by the electron. The dependence of these quantities is sensitive 
to their spatial variation inside the proton. 

Whereas the weak magnetism discussed above is a vector e-A^ interaction, an axial vector 
e-N coupling also exists, which is related to the proton's intrinsic spin. The parity violating 
e-p interaction depends on both of these quantities, and it is essential to determine the 
axial vector e-A^ form factor G\ in order to reliably extract G^. In general, G\ may be 
written as G\ = G\ + tiFa + -R^, where G\ is the contribution from a single Z-cxchangc as 
would be measured in neutrino-proton elastic scattering (Fig. lA), Fa is the nucleon anapole 
moment (1)^ and W is a radiative correction (Fig. IC). The constant r\ is ^-l^^Ow ^ '^'^^ 
where a is the fine structure constant. The anapole moment Fa is identified as the parity- 
violating coupling of a photon to the nucleon (Fig. IB) and is expected to be the largest of 
a class of higher order interactions, or radiative correct ions (^5/ It can arise as a result of, 
for example, a weak interaction between two quarks inside the nucleon. It is analogous to 
the nuclear anapole moment recently measured in atomic cesium (9)^ which is enhanced by 
parity violating interactions between nucleons in the cesium nucleus (10). Technically, the 
theoretical separation of the anapole moment from other radiative corrections W is gauge 
dependent, so the terminology associated with these amplitudes varies in the literature, 
and Fa itself is not cleanly defined theoretically. However, because the anapole moment is 
an electromagnetic interaction, it does not contribute to neutrino scattering, and is unique 
to parity- violating interactions with charged particles like electrons. Here we identify the 
observed difference as due to the anapole contribution (11). 

We performed a measurement of the parity violating asymmetry in the scattering of 
longitudinally polarized electrons from neutrons and protons (nucleons) in deuterium using 
the SAMPLE apparatus at the MIT/Bates Linear Accelerator Center. The apparatus was 
that used in our previously reported measurement on hydrogen (Ji.,5). In the deuterium 
experiment, the hydrogen target was replaced with deuterium, and borated polyethylene 
shielding was installed around the target to reduce background from low energy neutrons 
knocked out of the deuterium. Combining this measurement with the previously reported 
results allows an unambiguous determination of both the axial e-A'^ form factor G\ and the 
contribution of strange quarks to the proton's magnetic form factor 0%^. When presented in 
the context of the electron-quark couplings predicted by the standard model for electroweak 
interactions, our result also places new limits on the electron-quark axial couplings C-i^ and 

Parity violating electron scattering generally involves scattering of a longitudinally po- 
larized electron beam from an unpolarized target. The change in the count rate resulting 
from reversal of the beam polarization indicates a parity nonconscrving effect. For elastic 
electron scattering from a proton, the asymmetry in the count rate can be written as 

^ ^ / G^ \ eG\Gl + tG^^GI -e'(\-\ sin^ Q^) G^^G\ 

The three terms in Apy reflect the fact that it arises as a result of an interference between the 
electromagnetic and weak interactions, where j^f ordinary form factors associated 

with the proton's charge and magnetic moment. An equivalent expression can be written for 
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electron-neutron scattering. The kinematic factors r = £ = [1 + 2(1 + r) tan^(^e/2)] 

and e' = \J {}- — £^)r(l + r) can be adjusted to enhance the relative sensitivity of the exper- 
iment to the three contributions to Apv- 

In the SAMPLE experiment (12)^ the beam is generated from circularly polarized laser 
light incident on a bulk GaAs crystal, resulting in an electron beam with an average polar- 
ization of 36%. The electrons are then accelerated in durations of 25 /isec to an energy of 
200 MeV before encountering a liquid hydrogen (or deuterium) target. The polarization of 
the laser hght is reversed at the accelerator repetition rate of 600 Hz. After encountering 
the hydrogen target, the scattered electrons are detected in a large solid angle air Cerenkov 
detector covering angles between 130° and 170°. Figure 2 depicts 1 of 10 detector modules 
placed symmetrically about the beam axis. All backward scattered electrons with energies 
above 20 MeV travel faster than the speed of light in air and generate Cerenkov hght that 
is focussed onto a photomultiplier tube by an ellipsoidal mirror. The integrated photomul- 
tiplier current is proportional to the scattered electron rate, which is about 10^ s~^ during 
the beam pulse. A shutter in front of the photomultiplier tube allows measurement of back- 
ground coming from neutrons and charged particles. As described in ref. (4.), the pulse-height 
spectrum of the light-producing portion of the signal is characterized with dedicated runs 
using very low current beam. Prom these two pieces of information we determine the fraction 
of the signal due to the Cerenkov yield, which varies from mirror to mirror, but is typically 
60%. Several feedback systems are used to minimize the dependence of the properties of the 
electron beam, such as position, angle and energy, on its polarization state. 

The parity- violating asymmetry measured in SAMPLE is dominated by scattering from 
the neutrons and protons rather than the deuterium nucleus as a whole, and for our incident 
electron energy of 200 MeV it can be written as 

Ad = [-7.27 + 1.78G'^(T = 1) + 0.75^1^] ppm (3) 

(where "ppm" stands for "part per million"). The analogous expression for elastic electron 
scattering on the proton is 

Ap = [-5.72 + 1.55G\{T = 1) + 3.49(7^] ppm. (4) 

In these expressions we retain explicitly the isovector (T=l, or proton — neutron) compo- 
nent of G\. The small isoscalar component (proton -|- neutron) has been absorbed into 
the first constant term. The combined measurements allow an independent experimental 
determination of both quantities. 

The raw measured deuterium asymmetry is — 1.41±0.14 ppm, with a contribution due 
to helicity correlated beam parameters of 0.06±0.06 ppm. The measured background con- 
tribution to the raw asymmetry is 0.17±0.18 ppm, consistent with zero as expected. The 
dominant source of background is from photoproduced neutrons, which is known to have a 
parity- violating asymmetry of less than 1.5% of our measurement (13). In our analysis we 
assume the background asymmetry is negligible. Removal, on a detector-by-detector basis, 
of all dilution factors (mostly due to the beam polarization but also from background yields) 
and application of small corrections for contributing physics processes yields a value for the 
physics asymmetry of 



5 



-6.79 ± 0.64 (stat) ± 0.55 (sys) ppm . 



(5) 



The dominant sources of systematic error are the beam parameter corrections (5%), determi- 
nation of the beam polarization (4%), and determination of the signal-to-background ratio 
(4%). Uncertainties due to our lack of knowledge of the deuterium nuclear structure are 

about 3% (U ). 

Combining this measurement with the previously reported hydrogen asymmetry (5) and 
with the expressions in equations (3) and (4) leads to the two sets of diagonal bands in 
Fig. 3. The inner portion of each band corresponds to the statistical error, and the outer 
portion corresponds to statistical and systematic errors combined in quadrature. Our best 
experimental value of G%j is 

G'm{Q'^ = 0.1) = 0.14 ± 0.29 (stat) ± 0.31 (sys) . (6) 

In order to quote a result for the static contribution it is necessary to know the momentum 
dependence of G^. Wc use the calculation of (15), resulting in 

Us = [0.01 ± 0.29 (stat) ± 0.31 (sys) ± 0.07 (theor.)] /xjv (7) 

where we have added an additional third uncertainty (theor.) to account for error in the 
extrapolation coming from the range of the theoretical prediction. Combining this result 
with the known magnetic moments of the neutron and proton implies that the contribution 
of strange quarks to the proton's magnetic moment is —0.1 ± 5.1 %. 

The allowed region in Fig. 3 also provides a determination of the isovector axial e- 
N form factor. From a theoretical standpoint, the most uncertain contribution to G\ 
is the anapole term. The dominant contribution is expected to come from Fig. lA, 
and at the kinematic conditions of the SAMPLE experiment would result in a value of 
G\{T = 1)=— 1.071±0.005 (6,16). This value is analogous to that which would be mea- 
sured in neutrino scattering. Anapole effects have been estimated to substantially reduce 
the magnitude of G\{T — 1) (17), and a recent update of that theoretical treatment yields 
an expected value G\{T = 1) = —0.83 ±0.26 (18). This latter value is shown as the vertical 
band in Fig. 3. 

Our experimental result is 

G\(T = 1) = +0.22 ± 0.45 (stat) ± 0.39 (sys) (8) 

indicating that the substantial modifications of G\ predicted in (18) are not only present, 
but with an even larger magnitude than quoted. The deviation of our measured value of 
G\ from —1.071 can be interpreted as a large anapole moment of the nucleon. In the 
case of cesium, the anapole moment is enhanced by collective nuclear effects that are not 
present for a single nucleon. Understanding these higher order elcctrowcak processes is 
essential to reliably interpret precision studies of atomic parity violation and other tests of 
the electroweak theory. 

If further theoretical work fails to produce agreement with this experimental result, it 
may be necessary to consider the possibility that the discrepancy is due to new physics 
not contained in the standard model of electroweak interactions. In that regard, wc can 
recast our result in terms of the electron-quark axial coupling parameters C2u and C2d (6). 
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Present limits on these quantities are from the first parity violating electron scattering ex- 
periment (19) and from a measurement performed at the Mainz electron accelerator (20). 
The form factor G\{T = 1) is related to the combination C2u — C2d, resulting in our new 
experimental value of 



This represents a factor of three improvement over the precision of the Mainz measurement. 
In the Standard Model, the dominant contribution of Fig. lA would result in 



Adding in the corrections associated with the amplitudes displayed in Figs. IB and IC 

changes the value to —0.058 ± 0.02. Unless additional contributions to G\{T = 1) are 
identified that can improve the agreement, our measurement is in disagreement with the 
present standard electroweak prediction at the 1.5 cr level. 

An improved determination of will be provided with a new measurement of parity- 
violating quasielastic electron-deuteron scattering at a lower beam energy. The spatial de- 
pendence of the strange quark contributions and will be studied in detail in an upcoming 
program of parity violation measurements at the Thomas Jefferson National Accelerator 
Facility. 



C2u - C2d = +0.015 ± 0.032 (stat) ± 0.027 (sys) . 



(9) 




(10) 
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Figure 1 



FIG. 1. Feynman diagrams representing three contributions to the axial e-N couphng: A) sin- 
gle Z-exchange, B) parity violating photon exchange, which contributes to the nucleon's anapole 
moment, and C) a 7-Z box diagram typical of radiative corrections. 
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Figure 2 



FIG. 2. A schematic view of one module of the SAMPLE experimental apparatus. Ten mir- 
ror-phototube pairs are placed symmetrically about the beam axis. 
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Figure 3 

FIG. 3. A result of a combined analysis of the data from the previous two SAMPLE measure- 
ments. The two error bands from the hydrogen experiment (5) and the deuterium experiment 
are indicated. The inner hatched region includes the statistical error and the outer represents the 
systematic uncertainty added in quadrature. The ellipse represents the allowed region for both 
form factors at the la level. Also plotted is the estimate of the isovector axial e-A^ form factor 
G\{T = 1) obtained by using the anapole form ^tor and radiative corrections of Zhu, et al. (13). 



